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^ ■ ABSTRACT 

o ; 

Q^' We present millimeter-wave observations of several molecular ions in the disk around 

^ ■ the pre-main-sequence star DM Tau and use these to investigate the ionization fraction 
in different regions of the disk. New Submillimeter Array (SMA) observations of H2D+ 
J=li o - li 1, N2H+ J=4-3 and CO J=3-2 are presented. H2D+ and N2H+ are not 

\ detected and using the CO 3-2 disk size the observations result in an upper limit of <0.47 

\ K km s~^ for both lines, a factor of 2.5 below previous single-dish H2D"'" observations. 

' Assuming LTE, a disk midplane temperature of 10-20 K and estimates of the H2D"'" 

(-H ' o/p ratio, the observed limit corresponds to Njj2d+ < 4 — 21 x 10^^ cm~^. We adopt 

0^' a parametric model for the disk structure from the literature and use new IRAM 30 

Q . meter telescope observations of the H^^CO^ J=3-2 line and previously published SMA 

! observations of the N2H+ J=3-2, HCO+ J=3-2 and DCO+ J=3-2 hues to constrain 

, the ionization fraction, xi, in three temperature regions in the disk where theoretical 

\ considerations suggest different ions should dominate: (1) a warm, upper layer with 

I T>20 K where CO is in the gas-phase and HCO"*" is most abundant, where we estimate 

■ Xj ~ 4 X 10"^", (2) a cooler molecular layer with T = 16-20 K where N2H+ and DCO+ 

' abundances are predicted to peak, with ~ 3 x 10^^^, and (3) the cold, dense midplane 

. with T<16 K where Hjj' and its deuterated isotopologues are the main carriers of positive 

. charge, with < 3 x 10^^*^. While there are considerable uncertainties, these estimates 

\ are consistent with a decreasing ionization fraction into the deeper, colder, and denser 

\ disk layers. Stronger constraints on the ionization fraction in the disk midplane will 

\ require not only substantially more sensitive observations of the H2D"'" li^o ~ line, 

' but also robust determinations of the o/p ratio, observations of D2H+ and stronger 

\^ • constraints on where N2 is present in the gas phase. 

Subject headings: protoplanetary disks; astrochemistry; stars: formation; ISM: molecules; 
techniques: high angular resolution; radio lines: ISM 



1. Introduction 



The disks around pre-main sequence stars are the sites of planet formation. Models of disk 
evolution depend on many disk properties, including the total disk mass, the detailed density 
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and temperature structure, and accretion characteristics. Determining these properties are key to 
constrain how, when and where in disks planets can form. Since the main mass constituent in these 
disks is cold molecular hydrogen that cannot be observed directly, observations of dust and trace 
gas molecules are required to access the disk midplane properties. The ionization fraction Xi in disks 
is especially important because the magneto-rotational instability (M RI), believed to drive viscous 
accretion, requires ionization to couple the magnetic field to the gas (iGammid 119961 ). Insufficient 
ionization may suppress the MRI and create a quiescent "dead zone" with important implications 
for processes such as d ust settling, planetesimal formation, and planet migration (jCieslal 120071 : 
Matsumura et ahlbood ). 



According to staiidard chemical models, d isks are comprised of three characteristic layers (e.g. 



Willacy fc Millarjll998l : lAikawa fc Herbstll 19991 ). Figured] shows a schematic of the proposed layered 



structure focusing on ions. The uppermost layer is exposed to ultraviolet radiation and X-rays that 
dissociate molecules and photoionize molecules and atoms to produce, e.g. C"*". Deeper into the 
disk, the radiation field is attenuated and molecules such as CO becom e abundant. In this coole r 
molecular layer, HCO'*' is predicted to be the most abundant ion (e.g. lAikawa &: Nomural l2006l ) . 
At even lower heights, at temperatures below 20 K, CO begins to freeze out onto gr ains. N2 can 
remain in the gas-phase at slightly lower temperatures than CO (jOberg et al.ll2005l ). but in the 
deepest midplane layer, where temperatures drop below 16 K, the only molecules predicted to 
remain in the gas phase are H2, Hg" and their isotopologues. Indeed, II2D"'" and D2H"'" may provide 
the best, and possibly the on ly, observable tracers of the gas kinematics and the ionization fraction 
in the cold, dense midplane (jAsensio Ramos et al.l 120071 ). It is important to note that in addition 
to the vertical layers, each region has an inner truncation radius where the midplane temperature 
exceeds the CO and N2 freeze-out temperatures, respectively. Ions such as II2D"'" are therefore only 
expected in the outer d iisk, beyond th e CO snowline, whose radius will depend on the luminosity 
of the central star (e.g. lOi et al.l 120111 ). 



DM Tau is a pre-main sequence star located ~140 p c away from the Sun , surro unded by 
a gas-rich Keplerian circumstellar disk first discovered by iGuilloteau &: Dutreyl (119941). With a 



mass of 0.04 Mq, the disk is one of the more massive protoplanetary disks kno wn ([Andrews et al 



2011 



It is a trans ition disk with an imaged central hole of radius 19 AU (ICalvet et al.l l2005l : 



Andrews et al 



2011 



The DM Tau disk shows strong emission from species proposed to trace 



cold gas (jOberg et al.ll2010l ). in particular DCO"*" and N2 H"'", and has previo usly been detected 
in a variety of molecular lines using single-dish teles copes (IDutrev et al.lll996l ). It is therefore an 
especially good target to search for II2D"'" emission. ICeccarelli et al.l (|2004l ) reported a detection 
of the H2D+ li.o — li,i line at 372 GHz at the 4.7 a level t o ward DM Tau using the Caltech 
Submillimeter Observatory (CSO). However, IGuilloteau et al.l (j2006j) suggest that the DM Tau 
result is "at best a 2a detection" when analyzed with a proper model. To help to resolve these 
conflicting claims, we have used the Submillimeter Array to observe this H2D"'" line toward DM Tau 
with improved sensitivity. We analyze these results, together with observations of other molecular 
ions, to provide new constraints on the ionization fraction in different regions of the disk. 
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The paper is organized as follows. In §2, we describe the SMA observations and complementary 
IRAM 30 meter telescope observations of molecular ions in the DM Tau disk. In §3, we present 
upper limits on II2D+ and N2H+ J=4~3 emission and detections of the CO J=3-2 and H^^CO^ 
J=3-2 lines and other resolved observations of molecular ions from the literature. In §4, we use 
the line intensities together with an empirically derived disk structure to constrain the midplane 
temperature (§4.1) and to calculate upper limits on the o— Il2D^ column (§4.2), the total Il2D^ 
column (§4.3), the total midplane ion column (§4.4), the midplane ionization fraction (§4.5), and 
the ionization fraction in higher disk layers (§4.6). In §5, we compare the results with predictions 
from a generic disk chemistry model. 



2. Observations 

2.1. SMA observations of CO, N2H+ and 0-H2D+ 

DM Tau was observed with the Submillimeter Array on top of Mauna Kea, Hawaii, on 2010 
October 7. Six of the array antennas were available for these observations, arranged in a com- 
pact configuration spanning baselines of 16 to 77 meters. Using the dual receiver mode, we si- 
multaneously targeted the ortho-il2T)~^ li_o — li,i line at 372.4213 GHz and the N2H"'" J=4-3 at 
372.6725 GHz with the "high frequency" receiver, and the 345.796 GHz CO J=3-2 hue with the 
"low frequency" receiver. Table 1 summarizes the spectral line setup. At these frequencies, the 
synthesized beam sizes were 2'.'2 x 1'.'7 and 2'.'4 x 2'.'0, respectively. 

The observing conditions were generally very good, with atmospheric opacity r225 gHz ~ 0.04 
measured at the nearby Caltech Submillimeter Observatory, and stable atmospheric phase. This 
opacity corresponds to a zenith transmission of ~ 50% at the H2D"^ line frequency. Short observa- 
tions of the two calibrators 3C120 and J0449+113 were interleaved with the observations of DM Tau 
for gain calibration. The bandpass response was calibrated with observations of Uranus and the 
available bright quasars 3C273, 3C454.3 and 3C84. Observations of Uranus and Callisto provided 
the absolute scale for flux densities, via comparison to theoretical models for their emission, and 
the two calibrators agree within the expected uncertainties, but all reported data were calibrated 
with Callisto. The systematic uncertainty in the absolute flux scale is ~10%. The data were edited 
and calibrated with the IDL-based MIR software packag^. Continuum and spectral line images 
were generated and CLEANed using MIRIAD. To check that the flux calibration is accurate, we 
compared the co ntinuum fluxes of 0.22 Jy at 0.87 mm and 0.30 Jy at 0.81 mm with previously 



published values (I Andrews &: Williams! l2005l . 120071 ) and found that they agree within the reported 
uncertainties. 



http://www.cfa.harvard.edu/~cqi/mircook.html 
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2.2. IRAM 30m observations of U^^CO+ 

The J=3-2 line at 260.255 GHz was observed toward DM Tau on 2011 January 

25 with the IRAM 30 meter telescope. The observations were carried out with the EMIR 330 
GHz receiver, with a beam size of ~9'.'5 (FWHM). The receiver was connected to a unit of the 
autocorrelator with a spectral resolution of 320 kHz and a bandwidth of 240 MHz, equivalent to an 
unsmoothed velocity resolution of ~0.4 km s~^. Typical system temperatures were 300-400 K. The 
observations were carried out using wobbler switching with a 100" throw. Pointing was checked 
every ~2 hours on J0430+052 and J0316+413 with a typical accuracy of < 2". The main-beam 
brightness temperature was calculated from the antenna temperatures in pako using reported main 
beam and forward efficiencies (BeS a-nd Fgg ) of 88% and 53%, respectively. The data were reduced 
with the CLASS program, part of the GILDAS software pa^kagfl L.near ba^elmes were determined 
from velocity ranges without emission features, and then subtracted from the spectra. 



Results 



We did not detect either the H2D+ li^o — li,i line or the N2H"'" J=4-3 hue. Figure [2] shows 
imaging upper limits for these lines, along with moment maps for the CO J =3-2 line and thre e 
lines from other molecular ions, HCO+ J=3-2, DCO+ J=3-2 and N2H+ J=3-2 (loberg et ahlboiol ). 
The velocity gradient due to disk rotation appears similar in all of the moment maps, even for the 
weak N2H^ and DCO"'' emission. The apparent disk sizes are also similar within the uncertainties. 
To estimate the disk emission region, we fit an elliptical Gaussian to the CO J=3-2 visibilities. 
The resulting major and minor FWHM are 3'.'1 and 2'.'5. The effective disk size is defined as twice 
the FWHM of the CO emission, which will contain 98% of the emission assuming a Gaussian 
distribution. This provides a convenient definition of the disk size, which is not model dependent, 
but it is important to note that the defined disk radius does not directly correspond to the outer 
edge of the disk when the disk is modeled as a power law, or to the critical radius when the disk 
is modeled with the similarity solution. With the adopted definition, the CO 3-2 disk is 6'.'2x5'.'0, 
corresponding to a disk rad ius of ~400 AU, which compares well with previous gas size estimates 
(jGuilloteau &: Dutreylll994l ). The FWHM of all ions are consistent with the CO disk within the 
fit uncertainties. The synthesized beam sizes across the disk major axis vary between 2'.'0 and 
4'.'0, corresponding to spatial scales of 280-560 AU. The observations therefore only resolve the 
chemistry in the outer disk. 

Figure [3] shows the spatially integrated spectra of the species imaged in Figure [21 The inte- 
gration is carried out using an elliptical integration area the size of the CO disk. Table [2] lists the 
line fluxes integrated over the velocity range from 2 to 10 km s~^ The HCO+ J=3-2 has been pre- 
viously observed toward DM Tau using the IRAM 30m telescope by lDutrey et al.l (|l997l ). and the 



2http://www.iram.fr/IRAMFR/GILDAS 
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i ntegrated intensity of 4.1[0.5] Jy km s ^ compares well with the SMA value of 5.3[0.4] Jy km s ^ 



(|Qberg et al.ll20ld ). considering the mutual calibration uncertainties, which are in addition to the 



rms in brackets. 

To convert the SMA fluxes to intensities, we adopt the CO J=3-2 size scale for all of the 
emission lines and use the flux-to-intensity conversion formula: Tb = F x 13.6 x (A)^/(a x 6), where 
the constant factor is specific for the SMA observations, Tb is the intensity in K, F the fiux per beam 
(here the 6'.'2x5'.'0 CO disk) in Jy, A the wavelength of the line in millimeters and a and h are the 
major and minor axes of the beam (disk) in arc seconds (Table [2|). For the H2D"^ and N2H+ J=4-3 
lines, we calculate 2a upper limits with the standard formula: a = rms x FWHM/^Jn^, where 
rich is the number of channels across the CO 3-2 line at FWHM. The 2(T-limit on the H2D^ line 
intensity is <0.47 K km s~^. Averaged over the 22" beam of the Caltech Submillimeter Observatory, 
this 2a upper lirnit cor responds to <0.03 K km s~^, a factor of 2.5 below the claimed detection by 



Ceccarelli et al.l (j2004l ). This conclusively rules out the previous detection claim and entails that 
H2D"'" is not yet detected in any protoplanetary disk. This is consistent with the upper limits found 
by Chapillon et al. (subm. to A&A) using the JCMT. 

Figure [3] also includes the detection of H^'^CO^ J=3-2 from the IRAM 30 meter telescope. 
The reported intensity in Table[2]has been scaled to account for beam dilution by 9.5'^/(6.2 x 5.0), 
i.e. the ratio of the IRAM 30 meter beam area over the CO disk area. The fact that the H^'^CO^ 
line is detected indicates that the main isotopologue is optically thick, since if the emission ratio of 
HCO+/ H^3C0+ > 30, the 11^^00+ line would have been buried in the noise. 



4. The Disk Ionization Fraction 

We use the observations to constrain the ionization fraction in the disk in a four-step process. 
First, the observed upper limit on the o-H2D^ li^o-li,i line is used to constrain the 0-1120"*" column 
density N(o-H2D+) for a range of plausible temperatures for the disk midplane. Second, we use 
recent model results on the temperature dependent o/p ratio to convert N(o-H2D"'") limits to total 
{o + p) H2D+ limits, i.e. N(H2D"^). Third, we use calculations of the H;|"/H2D+/D2H+/D;|" ratios 
for the appropriate range of temperatures to constrain the total column density N(^H3_a;D2;) 
of these ions in the midplane. Finally, we use the total column density of ions in the midplane 
together with a physical model of the DM Tau disk and measurements of N(HCO"'"), N(DCO^), 
and N(N2H+) to constrain Xi in the disk midplane, an intermediate layer where CO has started to 
freeze out, and the upper warm molecular layer. 

Before embarking on this series of steps, we place observational constraints on the temperatures 
in the disk midplane with the aim to limit the range of temperatures to consider when deriving a 
midplane Xi from the 0-H2D"'" li,o-li,i- For this investigation, we define the midplane as the region 
where all heavy elements are frozen out and Hj^ and its isotopologues are the dominant charge 
carriers. Theoretically, the midplane temperature is limited by the freeze-out temperatures of 
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CO a nd N2, since both react readily with Hj^, resuhing in HCO"'' and N2H"'" ( Asensio Ramos et al 
2007^). Since N2 is maintained in the gas phase at sHghtly lower temperatures than CO (jOberg et al 



2005l ). the N2H+ rotational temperature provides the strongest observational limit available on the 



midplane temperature. 



4.1. Limits on the midplane temperature from N2H+ 10, 3 2 and 4 3 



To estimate the N2H"^ rotational temperature, we assume LTE and co-spatial emission from 
the N2H"'" 3-2 and 4-3 lines. This is a reasonable approximation because chemical models predict 



(Bergin et al. 


2006: 


Dutrev et al. 


2007) 



thin emission, unity filling factor within the CO 3-2 disk, and molecules thermalized at a single 
rotational temperature. 



v^^S J 



(1) 



from iThi et al.l (|2004l ). where N is the total number of molecules, Q{Trot) is the temperature 
dependent partition function, the energy of the upper level in K, and Trot the rotational tem- 
perature. We obtain the partition function Q(Trot) by second order interpolation of the listed 
values in CDM^, is gathered directly, and -Eu is calculated from the listed E'low values: 
Ea = -Eiow + 0.0334 X zv, where u is the line frequency in GHz. Then an upper limit on the 
N2H"'" temperature can be calculated from 



p4-3 l?3-2 

"-rot — — 7 — : T-r-r, — ^- 



In 



Using the integrated intensity of the J=3-2 line and the 2a upper limit on the J=4-3 line, we 
obtain a 2a upper limit on the N2H"'' rotational temperature of <26 K. This is consistent with the 
predictions that N2H"'' should become abundant below the CO freeze-out temperature of 20 K, but 
do es not provide any st ronger constraints. It is also consistent with the N2H"'' J=l-0 line reported 



by iDutrev et al.l ( 20071 ). We estimate the line flux from their Figure 1 to be <0.1 Jy km s~^. 
Using this limit and the integrated intensity of the J=3-2 line results in a lower limit on the N2H"'" 
rotational temperature of 11 K. More sensitive observations of the N2H"'" lines would be valuable 
to constrain this excitation temperature further. In the absence of more detailed information, we 
calculate LTE column densities of midplane ions for a range of temperatures, from 10 to 20 K. 



■^The Cologne Database for Molecular Spectroscopy at 



http://www.astro.uni-koeln.de/cdms/cataIog 
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The 20 K boundary is set by laboratory experiments on CO and N2 freeze-out - above this tem- 
pera ture a majority of CO and N2 are maintained in the gas-phase quickly destroying any formed 
Hj^ ( Bisschop et al.ll2006l ). The CO snowline was recently repo rted to correspond to a freeze-out 



temperature of 19 K, consistent with predictions (jOi et alj 12011 



4.2. N(o-H2D+) 

In general, column densities are calculated most accurately by constructing a self-consistent 
chemical-physical disk model and applying a radiative transfer code to predict line emission profiles 
that can be compared directly to observations. Without a detection of the H2D^ line, however, it 
is difficult to constr ain where the H2D"'" emis sion arises in the disk. We therefore opt for a classical 



LTE approach (e.g. Ivan Dishoeck et al.ll2003l ). calculated for the range of reasonable disk midplane 



temperatures (see ^4.ip to put limits on the disk averaged N(o-H2D"'"). 

A significant complication for the LTE calc ulation is that ortho jo) and para (p ) H2D"'" can 



behave as separate species at low temperatures (jFlower et al. 



2004 



Sipila et al.l 120101 ) resulting in 



IP) 



a non-thermal partitioning between the o and p states. To calculate the partition function for 
0-H2D+, we use a two-level approximation, which assumes no transitions between the o and p 
H2D~'', and that transitions between the ground state and first excited state dominate within each 
"species" of H2D"'". This approximation is expected to be valid at temperatures of 10 to 20 K, since 
Eu is 104 K for the first excited state of 0-H2D"'". The 0-H2D"'" partition function is calculated from 



Q = QuX e 



(3) 



where the degeneracies, Qu and gi, are equal to 3. 



Assuming LTE (Equation 1), excitation temperatures in the range 10 to 20 K and the partition 
function appropriate for these temperatures, N(o-H2D"'") limits are calculated from the observed 
upper limit on the H2D^" J=li^o-li,i line flux. As shown in Figure 3a, the limits on the disk average 
N(o-H2D"'") range from 1.3 
temperature). 



2.7 X 10 cm (with the highest value corresponding to the lowest 



4.3. Total (o + p) N(H2D+) 



The limit on the total N(H2D"'"), i.e. o + p, depends on the 0-H2D"'" column d ensity and the 
temperature dependent o/p ratio. This ratio has been modeled by Sipila et al. ( 201o l) from 4 to 20 K 
for a molecular hydrogen density of 10® cm~^ and an interstellar grain distribution. The derived o/p 
ratio is 0.5-0.1 for temperatures of 10-20 K, with the lowest ratio at 13-17 K. While the assumed 
grain size distribution and density are more appropriate for dense clouds than , disk midplanes, the 
grain size distribution mainly affects the temperature profile, and ISipila et al.l (|2010l ) find that the 
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H2D"'" o/p is nearly constant with increasing density. Thus we expect the values calculated with 
these assumptions should be valid for disk midplane conditions. 

Figure 3b shows the result of applying the literature o/p ratios to the previously calculated 
0-H2D"'" column density limits to obtain limits on N(H2D^) as a function of temperature. These 
limits range from 4 to 21 x 10^^ cm~'^. The presence of a peak value at 13 K is a direct consequence 
of the fact that the limit on N(o-H2D"'") decreases with temperature, while the o/p ratio has a 
minimum at 13-17 K. 



4.4. Total Midplane Ion column density: N(^ H3_a.D+) 



The limit on the total column density of ions in the midplane, 'H^.-x^t)-! depends on 

N(H2D^) and the ratio of N(H2D^) to N(^ H3_2,-D^). Like the o/p ratio, this ratio is expected 
to v ary with density a r id temperature, as well as other environmen tal properties such as grain 



size (iFlower et al 



2004 



Ceccarelli Dominik 



2005 



Sipila et al.ll2010l ). This ratio will also depend 
on whether or not CO and N2 are frozen out completely, or if low abundances of these specie s 



can be maintained in the gas-phase through non-thermal desorption (lAsensio Ramos et al.l 120071 ) . 
Caselli et al.l (j2008l ) modeled the ratios of all deuterated Hj^ isotopologues for molecular hydrogen 
density 10^ cm~^ and interstellar grains as a function of temperature. As mentioned prev iously, the 
disk midplane density is likely higher than assumed in these models. ISipila et al.l (|2010l ) found an 
increasing importance of Dg" with increasing density at 10 K, but it is unclear whether this effect 
is impo rtant also at higher t e mper atures. Another potential unknown is the abundance of H+, 
though ICeccarelli &: DominikI (120051 ) found that its contribution is negligible at the high densities 
of the disk midplane. 



With these caveats in mind, we use the ratios derived bv lCaselli et al.l (j2008l ) with maximum 
depletion, i.e. minimum amounts of CO and N2 in the gas-phase, to calculate upper limits on N(^ 
H3_2,.D+) in the disk midplane as a function of assumed temperature. Figure 3c shows the results, 
with limits that range from 1.5 to 8.5 x 10^^ cm~^, with the peak value at 13 K. The modification 
of the shape compared to the H2D^ curve is due to the fact that the N(H2D^)/N(^ H3_a;D^) 
ratio peaks near 13 K and then falls off quickly at higher temperatures. 



4.5. Ionization Fraction (xj) in the Midplane 



To derive limits on the midplane ion abun dances, we adopt the model of the DM Tau disk 
density and temperature structure described bv [Andrews et al.l (|201ll ) based on parametric fitting 
of the spectral energy distribution and resolved observations of millimeter-wave dust emission. This 
model does not include any chemistry, but rather provides an observationally constrained physical 
structure, which can be used to derive the gas masses in different temperature layers in the disk. 
This is needed to convert an ion column density into an ion abundance. We apply a standard 
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interstellar gas-to-dust ratio of 100 to convert the dust densities in this model to gas densities. To 
determine a midplane mass, we define the disk midplane as all material <16 K, the temperature 
for complete N2 freeze-out in most astro physical environments based on laboratory measurements 
( Qberg et al.l l2005l : iBisschop et alj 120061 ) . i.e. the region where it is most likely that Hjj" and its 
isotopologues become abundant. In this disk layer, the density weighted temperature is 11 K, and 
the disk averaged gas column density is 15 x 10^^ cm~^. Using this column density, and the 
upper limit on the ion column density in the midplane of 4 x 10^'^ cm~^ from Table 3 and Fig. 4, 
Xi < 2.7 X 10"^" uh. Alternatively, if the onset of CO freeze-out at 20 K defines the region where 
H2D"'" is abundant rather than N2 freeze-out at 16 K, then the density averaged ion temperature is 
13 K, and Xi < 4.5 x 10~^*^ uh. Table [3] summarizes the midplane (<16 K) ionization results and also 
shows how the ionization fractions compare if H2D^ is only assumed to be present in the midplane 
or in the entire T<20 K region. It is important to keep in mind that all of these calculations are 
based on model dependent values, and better limits on the midplane ionization abundance requires 
not only more sensitive o— H2D"'" observations, but also direct constraints on the o/p ratio and the 
relative abundances of different isotopologues. In particular the latter estimate may be off by an 
order of magnitude from different model predictions of the importance of Dj]" compared to 1120"*". 



4.6. Ionization in different disk layers 



It is interesting to compare the estimates of ionization fraction in different layers of the disk, 
since this may impact e.g. MRI activation and accretion flow. The mass of each layer, defined by 
a temperature interval, can be estimated from the same DM Tau physical model as was used in 
§4.5. Between 20 K, where CO begins to freez e out, and the CO photodissociation region, IICO+ 
is predicted to be the most abundant ion (e.g. IWillacvl 120071 ). Below 20 K, as CO starts to freeze 
out, N2II"'" becomes more abundant, since reactions with CO are the main destruction pathway of 
N2II"'". Below 16 K, CO and N2 are completely frozen out and Hg" and its isotopologues become 
the main charge carriers. 

The pCO '* ' abundance peaks betwe en the region where HCO^ dominates and the midplane 



(jWillacvl I2OO7I : lAikawa &: Nomural 120061 ) and therefore is assumed to coexist with N2II"'" in the 
intermediate layer. In this layer, DCO'^'/HCO"*' exceeds unity and the HCO^ contribution to the 
charge balance should b e small. This model result is supported by comparing DCO"*" line intensities: 
Guilloteau et al.l (|2006l ) detected DCO+ 2-1 at 0.29 ± 0.03 Jy km s'^ and the intensity ratio of the 
3-2/2-1 lines implies a rotational temperature of 11-23 K, assuming LTE. 

For the midplane and the low-temperature intermediate layer, the emission regions of the ions 
are expected to be well confined, and the assumption of LTE at a single rotational temperature 
should be a reasonable first approximation. For the warm molecular layer where HCO^ dominates, 
the assumption of LTE is not a priori a good approximation since this layer spans a temperature 
range from 20 to 100 K or higher. However, the derived HCO^ column density is not very sensitive 
within this range; it varies by less than a factor of three when using Eq. 1 over the full range of 
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plausible temperatures and the LTE approach retains utility. As discussed below, the outcome of 
this LTE calc ulation compares w ell with a radiative transfer calculation of HCO^ emission in the 
DM Tau disk (|Pietu et al.ll2007l l. 



We estimate the disk averaged N(HCO+), N(N2H+) and N(DCO+) using Equation 1, the 
observed intensities from Tabled and the mass weighted temperatures from the DM Tau structural 
model of 38 K for the CO dominated layer (>20 K) and 18 K for the intermediate layer where 
N(N2H+) and N(DCO'^) are most abundant (16-20 K). Table 3 lists the results and Figure [5] shows 
the resulting column densities together with the HjJ' limit (LTE at 11 K), within their respective 
temperature regimes (<16 K, 16-20 K, and >20 K), corresponding to the different disk layers. 
Note that N(HCO"^) is calculated from H^^CO^ assuming a standard ^'^C/^^C ratio of 70, since 
the high opacity of the main isotopologue underestimates N(HCO^) by a factor of seven. 

The calculation results in a disk averaged DCO^/HCO+ ratio of ^0.05. This is comparable 
to a previous estimate of 0.035 toward the disk around TW Hya (|Thi et al.ll2004l ). Both exceed 
the cosmic D/H ratio of ~10~^ by three orders of magnitude, demonstrating an efficient deuterium 
fractionation in t hese d isks. The derived ratio is an order of magnitude above the ratio found by 
Guilloteau et al.l (j2006l ). who used resolved HCO"'' 1-0 and single-dish DCO+ 2-1 emission and a 
radiative transfer calculation rather than LTE. The derived disk averaged N(HCO'*") are almost 
identical over the inner 400 AU of the disk: 7 x 10^^ cm~^ from their tabulated column density at 



300 AU and power law index, compared to 9 x 10 cm~ calculated here assuming LTE at 38 K. 
The different ratios must thus be due to different assumptions on the DCO"'' emission conditions. 
The HCO^ column density compares well with model predict ions of the out er disk (100-400 AU). 
It is for example consistent with the disk chemistry model bv IWillacvl (j2007l ) within a factor of 3. 



Dutrev et al.l (|2007l ) observed N2H+ 1-0 with the Plateau de Bure interferometer and derived 



N2H+/HCO''' of ~0. 02-0. 03. This ratio agrees well with the results in Table 3, which implies a 
disk-averaged ratio of 0.03. 

The disk averaged N(H2) within each temperature regime are 15, 2.4, and 2.6 x 10^^ cm~-^, 
and these values are used with the derived ion column densities to calculate disk averaged Xi in 
each layer. Figure [5] shows that the HCO'*' abundance and the midplane ion limits are comparable 
at ~ 3 — 4 X lO"^'^ njj, while the abundance of DCO"*" -|- N2H"'" is an order of magnitude lower at 
3 X 10 
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Overall, the derived ion abundances suggest a decreasing Xi with depth in the disk, which 
might be expected from attenuation of ionizing radiation. The assumption that DCO+ and N2H"'" 
alone are representative for the degree of ionization in the intermediate layer may be oversimplified, 
however. The 16 < T < 20 layer is a transition region where some HCO^ may still be present 
and Hj]" may already have started to become important. A combination of better limits on the 
midplane ions as well as more detailed modeling of the intermediate layer with DCO+ and N2H"'" 
will be needed to confirm the depth dependent trend on Xi. 
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Discussion 



In the previous section the disk averaged column densities and abundances were calculated 
based on observed line intensities and a density and temperature structure of DM Tau from the 
literature. This particular structure model has not been used as a basis for chemical modeling. In 
this section we therefore use predictions on the ion chemistry in generic disk models to compare 
our observational results with the current understanding of the chemical compositions in disks. 

Model predictions on the abui idances of H2D"'" and other ions iii protoplanetary disks have 



20071: iFogel et al. 



2011 



Aikawa & Nomura 


2006: 


Willacv 


2007: 



) for generic T Tauri disk structures. In particular, the models of IWillacv 
(|2007l ) include abundances and column densities for all of the ions considered in this study, and 
therefore we focus on comparing the observational results with these predictions. We examine 
specifically the disk averaged ionization fraction, the ionization fraction in different disk layers, and 
the relative abundances of different ions. 



Willacvl (|2007l ) based their chemical calculations on the disk structure model of lD'Alessio et al 
(j200ll ). representing a generic T Tauri disk. The assumed accretion rat e, disk mass etc ar e within 



200 



Most 



a factor of a few of the values that have been derived for DM Tau (jCalvet et al 
disk chemistry predictions are therefore expected to be valid for DM Tau, except for absolute 
ion column densities. The chemical modeliiig is b ased on the UMIST database, updated with 
multiple deuterated species (iRodgers &: Milladl 19961 ) . Within this model framework, IWillacyi (120071 ) 
investigated the importance of cosmic ray desorption and photodesorption on the disk chemical 
abundances, which has a large effect on the ion abundance ratios. 



An important caveat when comparing observations and model predictions, is that IWillacv 



(j2007l ) presents column densities as a function of radius with no simple procedure to calculate disk 
averaged column densities. Observations can therefore only be compared with model predictions at 
specific radii. Prom the DM Tau structure model the observed emission is dominated by the outer 
disk; beyond 50 AU the column densities of the different temperature layers are almost constant, 
which entails that the larger emission region of the outer disk will dominate. In particular more 
than 2 /3 of the ion emission is expected to or iginate outside of 250 AU if the ions are present across 
the CO disk. In the models bv IWillacvl (|2007l ) the ion column densities are constant within a factor 
of 2 at 100-400 AU, except for the relative importance of different isotopologues. It is therefore 
useful to compare model predictions in this reg ime with the ob served disk averaged values as long 
as this factor of two uncertainty is considered. IWillacvl (|2007l ) provides the most complete model 
results at 250 AU and this radius is used in all comparisons between model and observations. 
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5.1. Disk averaged ionization fraction 



The main sources of ionizing radiation in d isks are ultraviolet (UV) photons, cosmic rays, 
radioactive decay and X-rays (e.g. IWillacvl 120071 ) . and the ionization level therefore depends on 
a combination o f the inc i dent fluxes and the attenuation of external photons with dust and gas 
column density. IWillacyi (|2007l ) include both the interstellar radiation field and a stellar UV field 
that is 500 times stronger at 100 AU. Because of efficient absorption by dust and H2, the UV 
ionization is mainly important in the disk atmosphere and upper molecular layer. Cosmic rays are 
assumed to penetrate unperturbed to the disk midplane resulting in an ionization rate of 1.3 x 10~^^ 



Radioactive decay also contributes to the midplane ionization, at a rate of 0.6 x 10 



The effect of X-rays are not included in the model and the model may therefore underestimate rcj 
in lower disk layers. 



D+, HCO 



In this model, the total ion column density (H3 , H2D~'', D2H 
250 AU is 2 — 3 X 10^^ cm^^, corresponding to an average Xj of 2 — 3 x 10 
sistent with the observed abundances and limits. In the related disk model of Asensio Ramos et al 



DCO+, N2H+) at 
iiH. This value is con- 



(|2007l ). the ion abundance depends on the assumed cosmic ray flux, as expected. 



5.2. Vertical ionization gradients 



In disks, including the model from IWillacvl (|2007l ) the temperature increases both inward and 
upward in the disk, resulting in a vertically layered structure outside of an inner truncation radius. 
In the model, HCO"*" extends through the disk, while DCO"*" is only abundant outside of 50 AU and 
II2D"'' outside of 100 AU. As mentioned above, beyond 100 AU the different ion column densities 
are basically constant, justifying the comparison between disk averaged observations and model 
predictions at 250 AU. 

Because all UV photons and probably some of the X-rays and cosmic rays are absorbed in the 
outer disk layers, the ionization fraction should decrease betwe en the outer d isk layers and the disk 
midplane. At 250 AU, in the disk chemistry models (Fig. 3 in IWillacvl 120071 ) the ion abundance in 
the upper molecular layer, where CO is in the gas-phase and HCO"*" is the dominant ion, is a factor 
of 5 higher than in the midplane (~ 5 x 10^^^ compared to ~ 1 x lO^^" uh). In the intermediate 
region, where DCO^ and N2ll^ abundances peak, the ion abundance is about a factor of 2-3 below 
the CO layer ion abundance. This is qualitatively consistent with the results of the analysis of the 
DM Tau observations (see Figure 4). 



5.3. Specific ion column densities, abundances and ion ratios 



cm 



HCO+ and DCO+: The HCO+ column density predicted bv IWillacvl ([20071) is ~ 3 x 10^^ 
at 250 AU and is found to not depend on the different desorption mechanisms considered. 
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In contrast, the DCO"'' column density predictions vary by an order magnitude, resulting in a 
DCO^/HCO"*" ratio of 0.1-1. The higher DCO^/HCO^ ratios appear with desorption through 
cosmic ray heating, which maintains CO in the gas-phase in the disk midplane where the temper- 
atures are low enough for efficient deuterium fractionation and enhanced DCO"*" production. The 
observations indicate that the disk averaged DCO^/HCO"*" abundance ratio is below 0.1, which 
suggests that CO is confined to the warm >20 K molecular layer. The observational results are 
consistent with including UV photodesorption, which is efficient only in the warm molecular layer. 

N2H^ : The N2H"^/HC0^ abundance ratio is predicted to be in the range from 2 to 9 x 10^^, 
but the models with the highest values also predict DCO^/HCO^ ratios of unity and can be 
excluded. Compared to the value of 0.03 derived from the observations, the models are an order 
of magnitude too low. The reason for this discrepancy is unknown, demonstrating that even the 
simplest nitrogen-based chemistry remains poorly constrained. 

H2D~^ : N(H2D"'") is predicted to be 10^^ cm~^, with midplane abundance ~10^^^ uh. These 
values are well below the observational limits, and the reason is that Djj' is the dominant carrier of 
ch arge in the midplane in th ese models, 20-50 times more abundant than In the disk model 

bv lAikawa &: Nomural ( 2006l l. the Il2D^" column density is predicted to be significantly higher, 10^'^ 
cm~^, likely because multiple deuteration is not considered. Both models are thus consistent with 
the current limits, and better constraints on the chemical makeup of the midplane will require 
not only more sensitive observations of 1120"*", but also robust o/p determinations, and additional 
constraints on the abundance ratios of the isotopolgues. 



6. Conclusions 

1. A new SMA observation of the H2D+ li^o ~ li,i line toward the DM Tau disk results in a 
disk-averaged intensity upper limit of <0.47 K km s~^. 

2. Assuming LTE and a temperature dependent o/p conversion factor, the observation corre- 
sponds to an H2D"'" column density <4-21xl0^^ cm~^, for midplane temperatures of 10-20 K. 

3. Based on SMA and IRAM 30 meter telescope observations of additional molecular ions toward 
DM Tau and a previously published model of the disk structure, we estimate ionization 
fractions of Xj = 4 x lO"^'' in the molecular layer where T>20 K (based on H^^CO^), Xi > 
3 X 10^^^ in the 16-20 K layer where N2H^ and DCO^ are common, and < 3 x 10^^'^ in 
the midplane below 16 K where isotopologues dominate. The ion emission is expected to 
be dominated by the outer disk, beyond 250 AU, because of the larger emission area. 

4. The estimated ionization fractions depend on adopted model predictions of o/p ratios, deuter- 
ation levels and the dominating ions in different disk layers. Observational constraints on all 
three parameters together with deeper 1120"*" observations are therefore key to calculating 
more accurate ionization fractions in disks. 
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Fig. 1. — Schematic of the expected disk ion distribution as a function of disk layer temperature. 
The full lines represent the CO photodissociation front, the dashed lines the CO snowline, and the 
dotted line the transition region where the last heavy molecules freeze out. 
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Fig. 2. — Moment maps of the CO J=3-2 line and the molecular ion lines toward DM Tau observed 
with the SMA. The contour level (in Jy beam~^) are shown in the bottom right corner of each panel 
and correspond to 5o" for CO and HCO"*" and 2a for the other species. The color scale indicates 
the velocity field around the velocity centroid of 6.1 km/s. The synthesized beams are plotted in 
the lower left hand corners of each panel. 
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Fig. 3. — Spatially integrated spectra of CO J=3-2 and all searched for ion lines toward the DM 
Tau disk, including non-detections of the N2H"'" J=4-3 and H2D"^ li^o ~ li,i lines. The HCO"*" and 
CO lines have been scaled down for visibility. 
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Fig. 4. — Panel a) shows disk averaged N(o-H2D^) limits as a function of the assumed rotational 
temperature. Panel b) shows the limits on N(H2D+), where the temperature dependence is due 
to a combination of the calculated N(o-H2D+) limits and literature values on the temperature 
dependent o/p ratio. Panel c) shows the column density limits on N(^ H3_j;D+), which includes the 
temperature dependent relative abundances of the different isotopologues (as well as the previous 
temperature dependencies) . 
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Fig. 5. — Disk averaged ion column densities (upper panel) and ion abundances (lower panel) 
calculated at the density averaged temperature in the disk layers where Hg", DCO+ + N2H+ and 
HCO+ (here traced by H^^CO"*") are predicted in chemical models to carry most of the positive 

charge. The dashed line marks the transition between the midplane (where N2 freezes out) and 
N2H+ layer, and the dotted line marks the CO freeze-out temperature. The combined uncertainties 
from line fluxes and rotational temperature estimates arc generally smaller than the symbols, except 
for H^^CO"'' where the excitation temperature range is larger. 
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Table 1. SMA target spectral lines. 



Chunk 


Frequency range 


Channels 


Resolution 


Lines 


Frequency 


T 

up 




(GHz) 




(km s-i) 




(GHz) 


(K) 


USB:S08 


345.75-345.85 


256 


0.36 


CO 3-2 


345.7960 


33 


USB:S08 


372.37-372.48 


128 


0.66 


H2D+ ll,o-ll,l 


372.4213 


104 


USBiSll 


372.63-372.73 


128 


0.66 


N2H+ 4-3 


372.6725 


45 



Table 2. Line data with 2a upper limits and la uncertainties and continuum fluxes. 



Line 


/ Fdv 1 Jy km s'^ 


A / mm 


/7dw / K km s-i 


Fcota I naJy 


CO 3-2 


16.7[0.2] 


0.87 


5.43 [0.06] 


216[22] 


HCO+ 3-2 


5.3[0.4] 


1.12 


2.9[0.2] 


143[14] 


DCO+ 3-2 


0.7[0.1] 


1.39 


0.62[0.08] 


104[10] 


N2H+ 3-2 


1.0[0.1] 


1.07 


0.52 [0.06] 


143[14] 


N2H+ 4-3 


<1.8 


0.81 


<0.47 


303 [30] 


H2D+ ll,o-ll,l 


<1.8 


0.81 


<0.47 


303 [30] 


H.^^CO+ 3-2 


0.70(0.08] 


1.15 


0.22(0.03]'' 





* After correction for beam dilution. 



Table 3. LTE ion column densities and abundances in different disk layers. 





T Layer [K] 


T,ot [K] 


Nx [10^2 cm-2] 


nx [10-1° nn] 


0-H2D+ 


<16/<20 


ll'' / \t 


<2.4'' / <1.9'' 


<0.08'' / <0.06'' 


H2D+ 


<16/<20 


11 / 13 


<10 / <21 


<0.35 / <0.6 




<16/<20 


11 / 13 


<40 / <85 


<2.7 / <4.5 


N2H+ 


16-20 


18 


0.31±0.03 


0.14 ± 0.02 


DCO+ 


16-20 


18 


0.36±0.05 


0.16 ±0.02 


H( '0 


>2() 


:!8 


9. 1±().8 


:!..") ± ().:', 



Density weighted temperature in the midplanc layer when T<16 K. 
^ Density weighted temperature in the midplane layer when T<20 K. 



